Quenched dynamics in interacting one-dimensional systems: Appearance of current 
carrying steady states from initial domain wall density profiles. 
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We investigate dynamics arising after an interaction quench in the quantum sine- Gordon model 
for a one-dimensional system initially prepared in a spatially inhomogeneous domain wall state. 
We study the time-evolution of the density, current and equal time correlation functions using the 
truncated Wigner approximation (TWA) to which quantum corrections are added in order to set the 
limits on its validity. For weak to moderate strengths of the back-scattering interaction, the domain 
wall spreads out ballistically with the system within the light cone reaching a nonequilibrium steady- 
state characterized by a net current flow. A steady state current exists for a quench at the exactly 
solvable Luther-Emery point. The magnitude of the current decreases with increasing strength 
of the back-scattering interaction. The two-point correlation function of the variable canonically 
conjugate to the density reaches a spatially oscillating steady state at a wavelength inversely related 
to the current. 

PACS numbers: 71.10.Pm,37.10.Jk,05.70.Ln,75.10.Pq 



I. INTRODUCTION 

A fundamental question in the study of strongly cor- 
related systems concerns how a quantum many-particle 
system prepared in an initial state which is not an exact 
eigenstate of the Hamiltonian evolves in time, and under 
what conditions the system at long times thermalizes as 
opposed to reaching a novel athermal stated This ques- 
tion is particularly relevant now due to experiments in 
cold-atomic gases which provide practical realizations of 
almost ideal many-particle systems where the interaction 
between particles and the external potentials acting on 
them can be changed rapidly in time^ 

Motivated by this, there has been considerable the- 
oretical interest in studying the time-evolution of one- 
dimensional systems which are initially prepared in a spa- 
tially inhomogeneous state by the application of external 
confining potentials. Nonequilibrium time-evolution is 
triggered when the external potentials are rapidly turned 
off which may be accompanied with a rapid change in 
the interaction between particles. For example, the 
time-dependent density matrix renormalization group 
(TDMRG) has been used to study the time-evolution of 
a domain wall in the XXZ spin-chain^i the conformal 
field theory approach to study domain wall time evo- 
lution in the transverse-field Ising chain at the gapless 
pointy and the Algebraic Bethe Ansatz (ABA) to study 
Loschmidt echos for the XXZ chain for an initial domain 
wall state. 6 ABA has also been used to study geometric 
quenches i.e., the time-evolution arising after two spa- 
tially separated regions have been coupled together. 7 The 
dynamics of hard-core bosons after an initial confining 
potential was switched off was studied in Ref. [||. Here it 
was found that the initial energy of confinement resulted 
in the appearance of quasi-condensates at finite momen- 
tum. Time-evolution of an initial density inhomogeneity 
after an interaction quench at the Luther-Emery point 



was studied in Ref. ^ where a power-law amplification of 
the initial density profile was found. 

In this paper we study how a one-dimensional (ID) 
system prepared initially in a domain wall state corre- 
sponding to a density p(x — > ±oo) = ±po, (where x is 
the coordinate along the chain, and po is a constant) 
evolves in time after a sudden interaction and potential 
quench. The ID system is modeled using the quantum 
sine-Gordon (QSG) model which captures the low en- 
ergy physics of a variety of one-dimensional systems such 
as the spin-1/2 chain, interacting fermions with back- 
scattering interactions arising due to Umpklapp pro- 
cesses, and interacting bosons in an optical latticed 

The QSG model is integrable, its exact solution can be 
obtained using Bethe- Ansatz*^ While this property has 
been exploited to a great extent to understand equilib- 
rium properties of many ID systems, extending Bethe- 
Ansatz to study dynamics is a daunting task, especially 
for the time-evolution of two-point correlation functions. 
Thus there is a necessity to develop approximate meth- 
ods to study this model. 

Here we investigate the time-evolution of the QSG 
model semiclassically using the truncated Wigner ap- 
proximation (TWA) to which quantum corrections are 
added in order to set limits on its applicability. 12 More- 
over our parameter regime corresponds to an interacting 
bose gas whose density is initially in the form of a domain 
wall. We study how this initial state evolves in time as 
a result of a sudden switching on of an optical lattice, 
which may or may not be accompanied by an interaction 
quench. An optical lattice is a source of back-scattering 
interactions or Umpklapp processes, that tends to local- 
ize the bosons. Our aim is to understand how this physics 
affects the time-evolution of the domain wall state. Note 
that domain walls like the one we study here have been 
created experimentally by subjecting equal mixtures of 
S7 Rb atoms in two different hyperfine states to an exter- 
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nal magnetic field gradient^ Studying quantum dynam- 
ics in such systems may soon be experimentally feasible. 

One consequence of quenched dynamics in integrable 
models is that the system often does not thermalize, with 
the long time behavior depending non-trivially on the ini- 
tial state. Here we find that an initial state in the form of 
a domain wall evolves at long times into a current carry- 
ing state even in the presence of a back-scattering inter- 
action of moderate strength. Moreover, this net current 
flow has interesting consequences for the behavior of two- 
point correlation functions. The lack of decay of current 
found here is consistent with the fact that the dc con- 
ductivity of a ID system is infinite even in the presence 
of back-scattering or Umpklapp processes^ The origin 
of the infinite conductivity is the large number of con- 
served quantities in a ID system, where some of them 
have a nonzero overlap with the current ; i 15 thus pre- 
venting an initial current carrying state from decaying to 
zero. 

We also justify the steady state current obtained from 
TWA by studying the QSG model at the Luther-Emery 
point. The Luther-Emery point is an exactly solvable 
point in the gapped phase of the model. In particular 
we study how an initial current carrying state evolves 
with time and find that a steady state current (albeit of 
reduced magnitude) persists at long times. We also study 
how this current affects two-point correlation functions. 

Since the QSG model is a simplified model that ne- 
glects band-curvature and higher-order back-scattering 
or Umpklapp processes, an important question concerns 
to what extent it can capture quenched dynamics in real- 
istic systems. The nonequilibrium time-evolution of the 
above domain wall initial state was studied both for the 
exactly-solvable lattice model of the XX spin chain, and 
its continuum counterpart, the Luttinger models The 
study of the density and various two-point correlation 
functions revealed that both the lattice and the con- 
tinuum model reached the same nonequilibrium steady 
state, but differed in the details of the time-evolution. 
Continuum theories are far easier to handle both numer- 
ically and analytically than their lattice counterparts. 
Therefore to what extent they can capture the steady 
state behavior after a quench for general parameters is an 
open and important question which is beyond the scope 
of this paper. 

The paper is organized as follows. In section [TT] we 
study the time evolution of an initial domain wall state 
after a quench employing TWA. Results for the den- 
sity, current and two-point correlation functions are pre- 
sented. In section lTOl we present results for the first quan- 
tum corrections to TWA for some representative cases 
and discuss the general applicability of the TWA results. 
In section [IV] we present results for a quench at the ex- 
actly solvable Luther-Emery point for an initial current 
carrying state. Here results for the steady-state current 
as well as two-point correlation functions are presented. 
Section fVl contains our conclusions. 
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FIG. 1: (Color online) Density at time t=15 after the quench 
for 7f=0.9, 7=1 and several different g. 




FIG. 2: (Color online) Contour plots for the magnitude of 
the density for A"=0.9, 7=1 and for values of g a). g=0.05, 
b). g=0.2, c). g=0.6 and d). g=1.0. The density at t=0 is 
p(z)=(l/4) tanh (x/3). 




FIG. 3: (Color online) Contour plots for the -(current) 
for K=0.9, 7=1 and for values of g a). g=0.05, b). 
g=0.2, c). g=0.6 and d). g=1.0. The density at t=0 is 
p(as)=(l/4)tanh(a;/3) 
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FIG. 4: (Color online) The current at £=15 for g=0.05, 7=1 
and different K. 



FIG. 6: (Color online) Time evolution of the current after 
spatially averaging over a strip of width Sx = 40 centered at 
x = 0. 
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FIG. 5: (Color online) The current at £=40 for K=0.9,7=l 
and different g. 
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FIG. 7: (Color online) Dependence of the steady-state average 
current on interaction g for 7 = 1 and K = 0.9. 



II. TIME-EVOLUTION USING THE 
TRUNCATED WIGNER APPROXIMATION 

We start with an initial state which is the ground state 
of the Luttinger liquid, 
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where in terms of bosonic creation and annihilation op- 
erators bp, b p — j 
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and [4>(x), -d y 6(y)] = iS(x — y). Above, vp is the 
Fermi velocity or the velocity of the bosons, a a short- 
distance cut-off, p the momentum, L the length of the 



system, and h(x) is an external chemical-potential which 
couples to the density p(x)=—^d x (f>(x). In the ground 
state of Hi the density simply follows the external field 
(p(x)) = ^^;h(x). We choose h(x)=h tanh (x/£) so that 
the initial density is a domain wall of width £. We 
study the case where at time t=0 the external field h(x) 
is switched off. At the same time an optical-lattice is 
suddenly switched on which may be accompanied by an 
change in the interaction between bosons. Thus the time 
evolution for t > is due to the quantum sine- Gordon 
model, 



H f 



h 1 lh 



K {dj{x)f + -1 (d x cp{x)f 



g / dx cos (■~f(j>(x)) 



(4) 



Here u=vf/K, K being the Luttinger parameter and 
g the strength of the back-scattering interaction arising 
due to a periodic potential. The ground state of Hf has 
two well known phases^ the localized (gapped) phase 
characterized by ((f)) ^ 0, and a delocalized (gapless) 
phase. The periodic potential is a relevant parameter 
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for 2 — ^j— > 0, implying that the ground state has a 
gap for infinitesimally small g. On the other hand for 2 — 
->-£- < 0, a localized phase arises only for back-scattering 
strengths larger than a critical value (g > g c ). We will 
study quenched dynamics for parameters that are such 
that g is a relevant perturbation in equilibrium. Note 
that the initial domain wall state is not an exact eigen- 
state of Hf. Neither is it related to the classical solitonic 
solution of the QSG model since the latter is a domain 
wall in the <f) field,— while our initial state is a domain 
wall in d x (j). 

When 5=0, the time evolution of the system can be 
solved exactly^ For this case an initial density inhomo- 
geneity shows typical light-cone dynamics^ by spreading 
out ballistically in either direction with the velocity u, 
i.e., p(x, t)=2^~ [h(x + ut) + h(x — ut)]. Since the sys- 
tem is closed, tfie energy is conserved. However during 
the course of the time-evolution, the energy density is 
transferred from the density to the current, the latter 
having the form 

, x 1 96 
j(x,t) = -— (5) 

7T OX 

= ^no [h{ * x - ut) ~ hi * x+ut)] (6) 

In particular for h(x)=ho tanh(ai/£), the energy density 
at t=0 is £=ff (p(») 2 ~ hl/(2miK 3 ), while at long 
times, and for positions within the light-cone (ut > \x\ ) 
the energy density is £= 2±2 j^-j 2 where j=— ho/(iruK 2 ). 

Note that while any initial density profile will give rise 
to transient currents, the special feature of a domain wall 
density profile is that for a system of infinite length, the 
steady state behavior is characterized by a net current 
flow. In particular for any finite time the current flows 
across a length \x\ = ut of the wire connecting the regions 
of high and low densities ±p a — ±h /(Trvp) at the two 
ends. 

We now explore how the time-evolution of the density, 
and the long time behavior of the current and two-point 
correlation functions is influenced by a back-scattering 
interaction (g ^ 0). The results are obtained using TWA 
which involves solving the classical equations of motion 
with initial conditions weighted by the Wigner distribu- 
tion function of the initial state. Thus TWA is exact 
when g=0 while the effect of g is the leading correction in 
powers of hJ2. Since Hi is quadratic in the fields, it can be 
diagonalized by a simple shift, i.e., Hi—J2 p ^o v p\p\ a p £ h» 
where b p =a p + h v j (v p 2n\p\L) , h p being the Fourier 
transform of h(x). The initial Wigner distribution func- 
tion for the a p fields are Gaussian and are accessed by 
a Monte-Carlo sampling. This is followed by a Fourier 
transform defined in Eqns.[2]and[3]which gives the (j> and 
6 fields at the initial time t=0. The classical equations of 
motion are then solved on a lattice up to a time t. All the 
data sets presented here are accompanied with error bars 
associated with the Monte-Carlo averaging. Lengths will 
be measured in units of the lattice spacing a which is also 
set equal to the short-distance cut-off a. Energy scales 



will be in units of vp/a. The results will be presented for 
/i =7r/4 and an initial domain wall of width £=3. 



A. Time evolution of the density 

Fig. [1] shows the density at a time t=l5 after the 
quench for K—0.9, 7=1 and several different g. The do- 
main wall is found to broaden with time with a velocity 
which is reduced from the velocity of expansion u=vp/K 
when g = 0. Moreover, unlike purely ballistic motion, 
the shape of the domain wall changes during the time- 
evolution. The behavior of the density is clearer in the 
contour plots in Fig. [2] For small g, the time-evolution 
shows a light-cone behavior along with the appearance of 
spatial oscillations within the light-cone. The amplitude 
of the oscillations increase with g, while the wavelength 
of the oscillations is set by po. Increasing g gradually 
blurs the light-cone, and eventually for very large g the 
domain wall mass becomes so large that it hardly moves 
during the times calculated here. 



B. Time evolution of the current 

The current behaves in a manner complementary to 
the density and consistent with the continuity equation. 
Fig. [3] shows contour plots for the current for parame- 
ters that are identical to that for the density shown in 
Fig. [2] The current, like the density, fluctuates in space 
and time, but on an average reaches a non-zero steady 
state within the light cone for g values that are not too 
large. Fig. |4] shows the current at time £=15 for a given 
g and different K. As K decreases, the current increases 
as one expects from the analytic result for g=0. Fig. [5] 
shows how the current behaves for a fixed K and different 
g. Increasing g not only reduces the overall velocity of ex- 
pansion, but also reduces the magnitude of the current. 
Fig. [6] shows how the current spatially averaged over a 
strip of width 40 centered at the origin evolves in time. 
There is a clear appearance of a current carrying steady 
state whose magnitude decreases with g. Note that the 
spatial averaging under-estimates the time required to 
reach steady state as it under-estimates the amount of 
current for ut < 20. 

The dependence of the steady state current on g is plot- 
ted in Fig. [7] after time-averaging the current in Fig. [6] 
over a time window t = 40 — 100 in order to eliminate the 
temporal fluctuations. The steady state current is found 
to decrease linearly with g for g <C 1. Note that when 
g 7^ 0, the current does not commute with Hf. Yet the 
system reaches a current carrying steady state. This is 
due to the fact that the QSG model has a large number of 
other conserved quantities, some of which have a nonzero 
overlap with the current operator, thus preventing the 
current to decay to zero. The lack of decay of an initial 
current carrying state is also the origin of an infinite con- 
ductivity in many integrable systems*^ It was argued in 
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Ref. [I!] that at least two different non-commuting Ump- 
klapp processes are needed to violate conservation laws 
sufficiently so as to render the conductivity finite and 
thus cause the current to decay to zero. 

C. Steady-state correlation functions 

In this subsection we will study the following two equal 
time two-point correlation functions, 

C ee {xt;yt) = ( e «C*.*> e -*(».*)) (7) 
C H (xt;yt) = (e^O^e-**"'*)) (8) 

These are found to reach a nonequilibrium steady state 
for a time ut > \x\, \y\, i.e. for observation points that are 
within the light-cone. The result for Cee for K=0.9 and 
different g is plotted in Fig. [5] The TWA result for g=0 
is in agreement with the analytic result^ Cee(x,y;ut > 

\x\, |y|)=exp [ih (y - x)/(v F K)} (a/ \ x - y \f+ K ^/\ 
Thus when g=0, the correlation function decays as a 
power-law with a slightly larger exponent than in equi- 
librium (the latter being Cgg(x,y)—(a/ \ x — y \) 1 ^ 2K " 1 ). 
Moreover Cee shows oscillations at wavelength 
\= 2lTV h F g K =2 / j , j being the steady-state current 
within the light cone. Note that the results for ho = 
were obtained previously in Ref. Ojjl where the authors 
studied an interaction quench from a homogeneous ini- 
tial state. The physical reason for the spatial oscillations 
when ho ^ is the dephasing of the variable canonically 
conjugate to the density as the domain wall broadens. 
This implies a dephasing of transverse spin components 
in the XX spin chain resulting in a spin-wave pattern 
at wavelength Ap^& For a system of hard-core bosons, 
oscillations in Cee has the physical interpretation of 
the appearance of quasi-condensates at wave-vector 
k=2n/\£ 

The TWA results presented here show that these ef- 
fects can persist even in the presence of a back-scattering 
interaction, at least within the continuum model. In 
particular Fig. [5] shows that when g ^ 0, Cee retains 
the spatially oscillating form albeit at a wavelength that 
increases with increasing g. Just as for 5=0, one ex- 
pects the current to set the wavelength of the oscilla- 
tions. To check this Figs. |H] and [TU] show a comparison 
between C'ee(xt, yt) and Cee(ho=0)(xt, yt) cos(7rj(a; — y)) 
where Cee(ho=0)(xt 1 yt) is the correlation function at 
long times after a homogeneous quench from Hi to Hf 
(ho=0 in Hi), while j is the spatially averaged steady- 
state current in Fig. [7l The agreement is found to be 
good at least for small g. The second effect of g on Cgg 
is to give rise to a faster decay in position. This decay 
also becomes faster for a given g and on decreasing K 
(not shown here) which takes the system deeper into the 
gapped phase of the QSG model. 

Fig. HU shows the behavior of correlation func- 
tion after it has reached a steady state. The result for 
.9=0 is 16 ' 19 C^(x,y;ut>\xl\y\)=(a/\x-y\) (1+K2)/4 



which is also characterized by a slightly faster power- 
law decay than in equilibrium (the latter being 

^^>( x >y) = ( a / \ x — y \) K/2 )- Further unlike Cee, 
has no memory of the initial spatial inhomogeneity for 
3 = 0. However this is no longer the case for nonzero g. 
Fig. fTTI shows that for small g, C^ can also show spatial 
oscillations. Moreover, there is no appearance of long- 
range order until g is 0(1), where the appearance of the 
Ising gap corresponds to a nonzero asymptotic behavior 
of the two-point correlation function. It is also consistent 
that the appearance of the gap in C^ coincides with the 
value of g for which the domain wall is almost static in 

Fig. m 

In equilibrium, the cos j(f) interaction is a relevant per- 
turbation for 2 > j 2 K/AJ^ Thus the Hf parameters con- 
sidered here are those for which the ground state is the 
gapped Ising phase for g of any strength. Yet there is no 
signature of the gap in the quenched dynamics for g <C 1 . 
For this case the domain wall motion is ballistic, and the 
Cee correlations persist over longer distances than in the 
gapped Ising phase. Similar observations have also been 
made in the study of an interaction quench both in the 
bose-Hubbard models and for a system of interacting 
fermions^i where it was found that the system contin- 
ued to show light-cone dynamics and gapless behavior 
for parameters that correspond to the equilibrium gapped 
phase. 

The time-evolution of the density and current in the 
XXZ chain for an initial domain wall state was studied 
in Ref. [3| employing TDMRG. There it was found that 
while a current persists within the gapless phase, it de- 
cayed to zero in the gapped phase. This result is different 
from what we find here where the current persists in the 
gapped phase as long as g is not too large. There could 
be two reasons for this difference. Firstly the parameters 
7, g and K that we use here, do not correspond to the 
parameters of the XXZ chain. Secondly, it is possible 
that the irrelevant operators that are not retained in the 
continuum model modify the long-time behavior, even 
though this was not found to be the case at the exactly 
solvable XX point [J z = 0, K = 1, g = 0) J£ 



III. QUANTUM CORRECTIONS TO TWA 

An important question concerns the validity of TWA. 
It was shown in Ref. [TU that in writing the time-evolution 
of an interacting system as a Keldysh path integral, TWA 
is the leading correction in powers of K. One may there- 
fore check its validity by expanding the path integral in 
higher powers of H, and identify when these contributions 
become significant. We evaluate the first quantum cor- 
rection along the lines of Ref. HI]- Below we briefly outline 
the approach. 

The expectation value of an observable 0(x,p,i) to 
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FIG. 8: (Color online) The equal time Cee{0t;nt) correlation 
function at f=45 for K=0.9, 7=1 and different g. 
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FIG. 11: (Color online) The equal time C,^(0i; nt) correlation 
function at i=45 and for if=0.9, 7=1 and different g. 
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FIG. 9: (Color online) The equal time Cee{0t;nt) correlation 
function compared with correlation function for a homoge- 
neous quench (/in=0) and modulated by cos(njn) for t=45 
and g =0.03, A"=0.9, 7=1. 
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FIG. 10: (Color online) The equal time Cge(0t;nt) correla- 
tion function compared with correlation function for a ho- 
mogeneous quench (ho=0) and modulated by cos(7rjn) for 
t = 45,3=0.07, if =0.9, 7=1. 
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where WoCxojPo) is the initial Wigner distribution, and 
<V(x, p, t) = / dy (x - y/2| 6(x, p, t) |x + y/2) e*-*/* 
is the Weyl symbol of the operator O. In the QSG 
model, x — > 0(x), p — > II(a;) = ^d x 0(x), and J dx<ip — > 
j 2?0(x)2?n(a;). In Ref. [l2| the author implemented this 
correction by allowing a stochastic quantum jump in the 
momentum variable during the time evolution. This is 
done as follows: for each Monte Carlo step, we choose a 
set of initial conditions, weighted by the initial Wigner 
distribution, in accordance with TWA. For each set of 
initial conditions, wc select a random position, x n , and a 
random time, r. During the classical evolution, the field 
n n (£) = n(cc n ,i) is given a quantum kick at time r by 
shifting n„(r) ->■ U n (r) + £_{At) 1 / 3 , where £ is a ran- 
dom weight chosen from a Gaussian distribution of zero 
mean and unit variance, and At is a small time interval. 
Here, we take At equal to the integration time step size, 
At. This process of sampling t, x n , and £ is repeated 
for a given set of initial conditions. Thus the quantum 
correction to TWA is^ 



tNgj 3 



i[7^(r)](e 3 /3-£)0(</),n,i)), (10) 



where N is the number of spatial points. 

The results for the first quantum correction for 7=1 
and 7=2 are shown in Fig.[T2"land Fig . [T3l respectively. As 
expected, the larger the coefficient 7, the larger the quan- 
tum fluctuations in the </> field, causing TWA to break 
down sooner. We find TWA to work very well for 7=1 
up to times t = 15. On the other hand, for the same 
times, the quantum corrections for 7=2 are significant. 

It is also important to understand whether the steady- 
state current is a result of the truncation scheme. To 
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FIG. 12: (Color online) TWA and first quantum correction 
for the density (main panel) and equal time Cee correlation 
function (inset) for K=0.9, <?=0.05, 7=1 and t—15. 



FIG. 14: (Color online) TWA and first quantum correction 
to the current spatially averaged over the light-cone ut for 
if=0.9 and 7=1 and several different g. 




FIG. 13: (Color online) TWA and first quantum correction 
for the density (main panel) and equal time Cee correlation 
function (inset) for K=0.9, <?=0.05, 7=2 and f=15. 



check this we plot the current evaluated from TWA along 
with the first quantum correction in Fig.[14]for 7=1 and 
K = 0.9. The current is now spatially averaged over the 
non-interacting light-cone |x| < ut. The quantum cor- 
rection is found to enhance the current. This is expected 
on the grounds that TWA underestimates the quantum 
fluctuations, and therefore underestimates the extent of 
gapless behavior in the dynamics. 



IV. QUENCHED DYNAMICS AT THE 
LUTHER-EMERY POINT FOR AN INITIAL 
CURRENT CARRYING STATE 

The main result of TWA was that a current carrying 
state can persist even in the gapped phase of a model. 
In this section we will explore this physics at the Luther- 
Emery point of the QSG model. In particular we will 
study how an initial current carrying state evolves in 
time when the back-scattering interaction is suddenly 
switched on. We will also explore the long time behavior 



of two point correlation functions. 

The Luther-Emery point is a special point of the 
QSG model where the problem is rendered quadratic 
after refermionization in terms of left and right mov- 
ing fermions iI>l,r^ To see this, we rescale the fields 
in Eqn. ([I]) as <f>' = 7</>/2 and 9' = 28/j. Further if 
K = 4/7 2 , then Hf may be written as 



H'j = — iu I dx 



+ m I dx 



i>R(x)d x i/}R(x) - ipl(x)d x ^L(x) 



where m = gira and 

^r{x) = 



VR c -i\4>'(x)-e'M) 



V2 



^L_ e i[^'(x)+e'(x)] 



(11) 

(12) 
(13) 



VR,l are Klein factors to ensure the correct anticommu- 
tation relations among the fermions. 

We construct an initial current carrying state which is 
the ground state of the Hamiltonian 



dx 



i-)ipn(x) 
u 



-tpl(x)(d x 



i-)il) L (x) 
u 



(14) 



where 2/i is the chemical potential difference between 
right and left movers. We then study the time-evolution 
of this state for t > due to the Hamiltonian Hj (Eq. ITT]) . 
This Hamiltonian has a back-scattering interaction of 
strength m, and no applied chemical potential difference 
between right and left movers (/1 = 0). 

Defining ^a/lC^) = / l^e'^V'fi/iC^)) * ne initial state 
is characterized by the occupations 

(^ R {k)^ R (k)) = 6(-uk + n) (15) 
(V4»V>L(fc)> = 9(uk-ii) (16) 
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The current is denned as 



where 



j(x) = u ip R (x)ip R (x) - ij[(x)ip L (x) 



(17) 



Thus the initial state is characterized by a current density 

Jo = MA (18) 

Since the theory is quadratic, the time-evolution can be 
studied in terms of 

^(M) = ^*(*)/(M) + iM*0$(M) (19) 
iMM) = Mk)f*(k,t) + Mk)g(k,t) (20) 

where f(k,t) = cos(ojfet) — i sin(wfci) cos(20fc), t) = 
— i sin(wfc£) sin 2#/c, wfc=V«i 2 + w 2 fc 2 , tan(2^)=m/ (ufc). 

Using the above it is straightforward to work out the 
current at long times after the quench. The current 
reaches a steady state 



j = j Q - (m/x) tan 1 (j ir/m) 



(21) 



Note that this result is very similar to that obtained by 
TWA (Fig. [7]) and predicts that the steady state current 
decays linearly with m for m <C jo, while it decays as 
1/m 2 for large m. The persistence of an initial current 
carrying state even in t he g apped phase of a Hamiltonian 
was also found in Ref. [22|. Moreover in agreement with 
Ref. [22| we find that the steady state current in the limit 
of very small initial current jo <C m is found to scale as 
the cubic power of the initial current j cx j'q. 

We now turn to the evaluation of the steady-state gap 
and two point correlation functions. In terms of bosonic 
variables </>' = 7</>/2 and 9' = 26 /j, the gap is 



,2t^'(x,t)\ _ 



(22) 



while the basic two-point correlation functions are 



= (^ R (xt)M^[(0t)i; R (0t)) (23) 

Cg,g,(x,t) = ( e -2^'( a =,t) e 2ifl'(0,i) ) 

= (^ R (xt)^ L (xt)ip L (0t)M0t)) (24) 
For long times after the quench we find 

u 2 /a 2 



,2i<j>'(x,t)\ 



-in 



In 



m 2 + it 2 jo 



A{a,m,jo) (25) 



where a is a short-distance cut-off. Thus the steady-state 
gap depends on the initial current jo- 

The two point correlations at long times are 

Cppfat) = \A(a, m, j )\ 2 + \^6(x)+ zl b + I a \ 2 (26) 

CW(x,t) = (^S(x) + iI b + I a j (-5(x)-iI b -I a 



iL 



(27) 



h 
h 

Ic 

Id 



» /U dk „ , U 2 k 2 

cos(fca;)- 



00 dk 



u 



2tt 



ka sin kx 



2 k 2 

21,2 



u A k 



^ u dk 
27 



sin(fea;) 



m 2 + u 2 k 2 
muk 



u 



2 k 2 



dk 



muk 



fi/u 



e ka ^- coskx— — - 

Z7r m z + u z k z 



(28) 
(29) 

(30) 
(31) 



For jo = (ji = 0), reduces to the expression 

derived in<^ For jo 7^ and long distances fix/u 3> 
1, mx/u > 1 we find, 



C^,f,'{x,t) = \A(a,m, j )\ 2 + 



1 



(32) 



Thus the correlations are found to decay very slowly (as 
l/x 2 ) in position to their long distance value of the square 
of the gap. This should be contrasted with the equilib- 
rium result where the decay to the long distance value 
is exponential.— It is also interesting to compare this re- 
sult with that of an interaction quench from an initial 
state which is the ground state of H-(fi = 0)22- For this 
case the decay to the long distance value is a power law 
(l/cc 6 ), but with a larger exponent than found here for 
the current carrying state. 

The expression for Ce>e' at long times after the quench 

is 



1 



Ce>e>(x, Moo) = — 



- 2 J 2 



m 2 + it 2 j% 



^-2iirj x/u 



(33) 



and shows a similar slow decay as 1/a; 2 in position (in 
contrast to an exponential decay in equilibrium). More- 
over, the current flow imposes spatial oscillations at a 
wavelength which is determined by the current. 

This quench at the Luther Emery point did not involve 
a change in the Luttinger parameter K. Significantly dif- 
ferent physics can occur after a similar quench that also 
changes the value of K. In Ref. the authors showed 
that changing K can lead to the existence of "super soli- 
tons" at the Luther Emery point, where initial density 
inhomogeneities spread out with amplitudes that grow 
in time. 



V. CONCLUSIONS 

In summary, we have performed a detailed study of 
quenched dynamics in an interacting ID system prepared 
initially in a domain wall state. The model, being in- 
tegrable, never thermalizes with the system reaching a 
nonequilibrium current carrying state which is robust 
even in the presence of moderate back-scattering inter- 
actions. The current has interesting consequences for 
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the correlation functions, most notably the appearance of 
spatial oscillations in the Cgg correlation function. Our 
predictions for the current can be tested experimentally 
using presently available one-dimensional optical lattice 
techniques'^ 
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